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O . We study the stabihty of cosmological scahng solutions within the class of spa- 

tially homogeneous cosmological models with a perfect fluid subject to the equation 
of state p-y = — 1)^7 (where 7 is a constant satisfying < 7 < 2) and a scalar 
field with an exponential potential. The scaling solutions, which are spatially flat 
isotropic models in which the scalar fleld energy density tracks that of the perfect 
fluid, are of physical interest. For example, in these models a signiflcant fraction 
5^ I of the current energy density of the Universe may be contained in the scalar field 

whose dynamical effects mimic cold dark matter. It is known that the scaling so- 
lutions are late-time attractors (i.e., stable) in the subclass of flat isotropic models. 
We flnd that the scaling solutions are stable (to shear and curvature perturbations) 
in generic anisotropic Bianchi models when 7 < 2/3. However, when 7 > 2/3, and 
particularly for realistic matter with 7 > 1, the scaling solutions are unstable; es- 
sentially they are unstable to curvature perturbations, although they are stable to 
shear perturbations. We briefly discuss the physical consequences of these results. 

1. Introduction. 

Scalar fleld cosmological models are of great importance in the study of the early 
Universe. Models with a variety of self-interaction potentials have been studied, 
and one potential that is commonly investigated and which arises in a number of 
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physical situations has an exponential dependence on the scalar field [1 - 5] . There 
have been a number of studies of spatially homogeneous scalar field cosmological 
models with an exponential potential, with particular emphasis on the possible 
existence of inflation in such models [1]. 

These models may also be important even if the exponential potential is too steep 
to drive inflation. For example, there exist 'scaling solutions' in which the scalar 
held energy density tracks that of the perfect fluid (so that at late times neither 
fleld is neglible) [2] . In particular, in [3] a phase- plane analysis of the spatially flat 
Friedmann-Robertson- Walker (FRW) models showed that these scaling solutions are 
the unique late-time attractors whenever they exist. The cosmological consequences 
of these scaling models have been further studied in [4]. For example, in such 
models a signiflcant fraction of the current energy density of the Universe may be 
contained in the homogeneous scalar field whose dynamical effects mimic cold dark 
matter; the tightest constraint on these cosmological models comes from primordial 
nucleosynthesis bounds on any such relic density [2 - 4] . 

Clearly these scaling models are of potential cosmological significance. It is conse- 
quently of prime importance to determine the genericity of such models by studying 
their stability in the context of more general spatially homogeneous models. It is 
this question that we shall address in this paper. 

II. The Scaling Solution. 

The governing equations for a scalar field with an exponential potential 



where Vq and k are positive constants, evolving in a flat FRW model containing a 
separately conserved perfect which satisfles the barotropic equation of state 



where the constant 7 satisfles < 7 < 2 (although we shall only be interested in 
the range < 7 < 2 here), are given by 



V = Voe 



-K(t> 



(1) 



P7 = (7 - 1)P- 



'7' 



(2) 
(3) 
(4) 



4) = -3H4> + kV, 



subject to the Friedmann constraint 



H' = lip, + l^' + V), 



(5) 
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where H is the Hubble parameter, an overdot denotes ordinary differentiation with 
respect to time t, and units have been chosen so that SttG = 1. We note that the 
total energy density of the scalar field is given by 

P4> = + V. (6) 



Defining 



_ 
X = 



y = 



V6H ' VSH' 
and the new logarithmic time variable r by 

dr 



(7) 



dt 



H, 



(8) 



equations (2) - (4) can be written as the plane-autonomous system [3] : 
x' = -Sx + ^^ny'' + ^x[2x'' + 7(1 - x^ - y% 



y =-^y 



-\l -Kx + 2x^ + 7(1 -x^ - y^) 



(9) 
(10) 



where a prime denotes differentiation with respect to r, and equation (5) becomes 

1^ + ^0 = 1, (11) 



where 



Pi 



which implies that < x^ + y^ < 1 for > so that the phase-space is bounded. 

A qualitative analysis of this plane-autonomous system is given in [3] . The well- 
known power-law inflationary solution for < 2 [1] corresponds to the equilibrium 
point X = Av/\/6, y = (1 - n'^/ey/^ {Q^ = 1, = 0) of the system (9)/(10), which 
is shown to be stable (i.e., attracting) for < 37 in the presence of a barotropic 
fluid. Previous analysis has shown that when < 2 this power-law inflationary 
solution is a global attractor in spatially homogeneous models in the absence of a 
perfect fluid (except for a subclass of Bianchi type IX models which recoUapse) . 

In addition, for 7 > there exists a scaling solution corresponding to the equi- 
librium point 



X = Xq 



^Jll, y = yo = [3(2-7)7/2«']^ 



(13) 
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whenever > 87. The hnearization of system (9) / (10) about the equihbrium point 
(13) yields the two eigenvalues with negative real parts 

(2 - 7) ± ^ V(2 - 7)[247^ - «^(97 - 2)] (14) 

when 7 < 2. The equilibrium point is consequently stable (a spiral for > 
247^/(97 — 2), else a node) so that the corresponding cosmological solution is a 
late-time attractor in the class of fiat FRW models in which neither the scalar-field 
nor the perfect fluid dominates the evolution. The effective equation of state for 
the scalar field is given by 

^ {p<f>+P^) _ 2a;g _ 
Pcl> + Vo 

which is the same as the equation of state parameter for the perfect fluid. The 
solution is referred to as a scaling solution since the energy density of the scalar field 
remains proportional to that of the barotropic perfect fluid according to fl/fl^ — 
«;^/37 — 1 [2]. Since the scaling solution corresponds to an equilibrium point of the 
system (9)/(10) we note that it is a self-similar cosmological model [6]. 

III. Stability of the Scaling Solution. 

Let us study the stability of the scaling solution with respect to anisotropic and 
curvature perturbations within the class of spatially homogeneous models. 

A. Bianchi I models. 

In order to study the stability of the scaling solution with respect to shear per- 
turbations we shall flrst investigate the class of anisotropic Bianchi I models, which 
are the simplest spatially homogeneous generalizations of the flat FRW models 
which have non-zero shear but zero three- curvature. The governing equations in 
the Bianchi I models are equations (3) and (4), and equation (5) becomes 

H' = l{p. + l<i>' + V^+^', (15) 

where = ^Egi?"^ is the contribution due to the shear, where Sq is a constant 
and R is the scale factor. Equation (2) is replaced by the time derivative of equation 
(15). 

Using the definitions (7), (8) and (12) we can deduce the governing ordinary 
differential equations. Due to the term in (15) we can no longer use this equation 
to substitute for in the remaining equations, and we consequently obtain the 
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three-dimensional autonomous system: 



x' = -Zx + ^ 3 + (7 - 2)n - 2y\ (16) 



y' =\y^-^\KX + 2+{^-2)^- 2/ ) , (17) 
1]' = 31){(7-2)(1)- 1) -2y2}, (18) 
where equation (15) yields 

1 -O-x^ -1/2 = s^iy-^ > 0, (19) 
so that we again have a bounded phase-space. 

The scaling solution, corresponding to the flat FRW solution, is now represented 
by the equilibrium point 

87 

x = xq, y = yQ, Q = l j. (20) 

The linearization of system (16) - (18) about the equilibrium point (20) yields three 
eigenvalues, two of which are given by (14) and the third has the value —3(2 — 7), 
all with negative real parts when 7 < 2. Consequently the scaling solution is stable 
to Bianchi type I shear perturbations. 

B. Curved FRW models. 

In order to study the stability of the scaling solution with respect to curvature 
perturbations we shall first study the class of FRW models which have curvature 
but no shear. Again equations (3) and (4) are valid, but in this case equation (5) 
becomes 

H^ = l{p^ + l^^ + V) + K, (21) 

where K = —kR~'^ and A; is a constant that can be scaled to 0, ±1. Equation (2) 
is again replaced by the time derivative of equation (21). 

As in the previous case we cannot use equation (21) to replace p^, and using the 
definitions (7), (8) and (12) we obtain the three-dimensional autonomous system: 



7-l)n+l{l + 2x^-y') , (22) 



/ o /3 2 3 
X = -3x + d -Ky + -X 

{ -^^x + (^7 - ^ j O + ^(1 + 2x^ -y')}, (23) 
^-l){Q-l) + l{2x'-y')\, (24) 
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where 

1-n-x'^ -y^ = KH-^. (25) 

The phase-space is bounded for k = or k = —1, but not for k = +1. 

The scahng solution again corresponds to the equihbrium point (20). The hn- 
earization of system (22) - (24) about this equihbrium point yields the two eigen- 
values with negative real parts given by (14) and the eigenvalue (87 — 2). Hence 
the scaling solution is only stable for 7 < | . For 7 > | the equilibrium point (20) 
is a saddle with a two-dimensional stable manifold and a one-dimensional unstable 
manifold. 

Consequently the scaling solution is unstable to curvature perturbations in the 
case of realistic matter (7 > 1); i.e., the scaling solution is no longer a late-time 
attractor in this case. However, the scaling solution does correspond to an equilib- 
rium point of the governing autonomous system of ordinary differential equations 
and hence there are cosmological models that can spend an arbitrarily long time 
'close' to this solution. Moreover, since the curvature of the Universe is presently 
constrained to be small by cosmological observations, it is possible that the scaling 
solution could be important in the description of our actual Universe. That is, not 
enough time has yet elapsed for the curvature instability to have effected an appre- 
ciable deviation from the flat FRW model (as in the case of the standard perfect 
fluid FRW model). 

Hence the scaling solution may still be of physical interest. To further study its 
significance it is important to determine its stability in a general class of spatially 
homogeneous models. We shall therefore study the stability of the scaling solution 
in the (general) class of Bianchi type Yllh models, which are perhaps the most 
physically relevant models since they can be regarded as generalizations of the open 
(negative-curvature) FRW models. 

C. Bianchi VII/j models. 

The Bianchi VH/j models are sufficiently complicated that a simple coordinate 
approach (similar to that given above) is not desirable. To study Bianchi VH/j 
spatially homogeneous models with a minimally coupled scalar field with an expo- 
nential potential and a barotropic perfect fiuid it is best to employ a group-invariant 
orthonormal frame approach with expansion-normalized state variables governed by 
a set of dimensionless evolution equations (constituting a 'reduced' dynamical sys- 
tem) with respect to a dimensionless time subject to a non-linear constraint [6], 
generalizing previous work in which there is no scalar field [7] and in which there is 
no matter [8]. 

The reduced dynamical system is seven-dimensional (subject to a constraint) [9]. 
The scaling solution is again an equilibrium point of this seven-dimensional system. 
This equilibrium point, which only exists for > 87, has two eigenvalues given 
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by (14) which have negative real parts for 7 < 2, two eigenvalues (corresponding 
to the shear modes) proportional to (7 — 2) which are also negative for 7 < 2, 
and two eigenvalues (essentially corresponding to curvature modes) proportional to 
(37 — 2) which are negative for 7 < | and positive for 7 > | [9]. The remaining 
eigenvalue (which also corresponds to a curvature mode) is equal to 87 — 4. Hence 
for 7 < I {k? > 37) the scaling solution is again stable. However, for realistic 
matter (7 > 1) the corresponding equilibrium point is a saddle with a (lower) four- 
or five-dimensional stable manifold (depending upon whether 7 > 4/3 or 7 < 4/3, 
respectively) . 

IV. Discussion. 

Perhaps these stability results can be understood heuristically as follows. From 
the conservation law the barotropic matter redshifts as R~^^ . In subsection III. A 
we saw that in this case the shear redshifts as and so always redshifts faster 
than the matter, resulting in the stability of the scaling solution. We note that the 
bifurcation that occurs at 7 = 2/3 in subsection III.B corresponds to the case in 
which the curvature K is formally equivalent to a barotropic fluid with 7 = 2/3, 
and in which both the matter and the curvature redshift as R~^. For 7 > 2/3, the 
barotropic matter redshifts faster than and the curvature eventually dominates. 
A complete qualitative analysis of cosmological models with a perfect fluid and a 
scalar fleld with an exponential potential will be undertaken in future work [cf. 9]. 
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